Theoretical calculation of the surface energy of water 
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The estimation of the surface tension of water is theoretically dealt on the basis of the dipole 
molecular model. It is known that the experimentally determined surface tension of freshly ex- 
posed surface has a higher value than the nominal value of 73 [mN/m]. We calculated the value 
corresponding to the fresh surface where the reorientation of the molecules has not occurred. 

PACS numbers: 68.37.-d 



I. INTRODUCTION 

The surface tension of water, 7 = 73 [mN/m], is con- 
sidered to be a static value after the equilibrium of the re- 
arrangement of the molecules at surface is attained. The 
structure of surface orientation |lj of water molecules is 
not yet well-defined, though the simulations using molec- 
ular dynamics 0] have been studied. This may be one 
reason that the theoretical calculation of surface tension 
of water which explains the static value has not yet been 
reported. Theoretical treatments of the surface tension 
have been done for nonpolar van der Waals molecules by 
statistical mechanics, or by using intermolecular forces 
involving Hamaker constant on the basis of the Lifshitz 
theory [f, 0. 

These results showed considerably good agreement 
with experimental values of nonpolar materials. How- 
ever, theoretical calculation for the polar molecules seems 
to be hampered by the difficulty of evaluation of surface 
orientation. On the other hand, the surface tension of 
water is known to have higher value than the static one 
for freshly created surface. This has been shown by the 
dynamical measurement by the oscillating jet method 
where a value of 180[mN / m] is observed immediately af- 
ter the surface is created, and gradually lowers to its 
static value within less than 10 [mS] 6]. 

The higher value is thought to originate from the 
freshly exposed water surface where the molecules are 
not completely reoriented. We here focused on the cal- 
culation of the unoriented surface tension of water based 
on the dipole interaction model. The calculation of the 
contribution of additional terms of intermolecular forces 
are also presented. 

The surface energy which is defined as the (positive) 
energy per area of the surface, is equivalent to the en- 
ergy for creation of the surface. The surface energy orig- 
inates from the binding energy between molecules. One 
molecule inside bulk is surrounded by other molecules. 
Let us define the number of such molecules (number of 
nearest neighbour) as Ni, and the mean binding energy 
W (< 0) between nearest two molecules. Then the bind- 



ing energy per one molecule is 

Bk = -fW. 

On the other hand, the molecule on the surface is sur- 
rounded by less neighbours N s . Then the binding energy 
per molecule is 



-W. 



Therefore the molecule on surface has larger energy 
than that inside bulk. 



AE 



N s - N,; 



-W > 0. 



By using N[m 2 ] as molecule number density on sur- 
face, we define the surface energy per area. 



N s -Nj „ rAT 
u = WN. 



In the above discussion we considered only the nearest 
neighbour interaction. However, we must consider the 
long range force like Van der Waals and dipole-dipole 
interactions. For this purpose, our method for calculating 
the surface energy is as follows. 

1. Consider the virtual flat surface in bulk which di- 
vide bulk into two pieces. 

2. Calculate the binding energy of two molecules. 

3. Sum up the binding energies for all molecule-pairs 
across the virtual surface in bulk. 

4. The binding energy obtained above divided by two 
times surface area gives the surface energy per area. 

We will calculate the surface energy of water due to 
above procedure in the following sections. 

The surface tension 7 which we observe is defined by 
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7 = It — Ts, 



(1) 



2 



where u is the surface energy (per area), and s is the 
entropy density per area. This shows our calculation does 
not directly connect to observable 7. However, the sur- 
face entropy term is far possible to take into account our 
consideration, since it deeply relates with surface orien- 
tation. We just assume here the surface entropy term is 
much smaller than the surface energy u. Then the surface 
tension is approximated by surface energy and we com- 
pare the experimental value of surface tension of water 
to the surface energy which we calculate. 

Before starting our program, let us consider the dimen- 
sional analysis. The molecule of water has its own dipolc 
moment. |7J 



fj, = 6.471 x W- M [Cm] 



The molecular density is 



n = 3.35 x 10 28 kr 3 l 



(2) 



(3) 



By using these two quantities: /x and n, we obtain only 
one quantity which has dimension of energy density. 



2 

— — x n 16 



130[mJ/m 2 



(4) 
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where we have utilised the value -r^— = 9.01 x 10 9 . The 

47TC0 

surface tension of water is the same order as this value. 
This shows such a consideration has possibility to express 
the nature of water. 

II. DIPOLE-DIPOLE INTERACTION 

The Interaction energy between two electric dipoles 
has the form. 



v = 1 ■ $b 3 (Pa ■ r)(pB ■ r) , ^ 
47reo r 3 r 5 



Fig. 1 

By using the coordinates 9 a, Ob, 4>a, 4>b represented 
in Fig.l, This energy is written in the form. 



V 



[sin 9a sin Ob cos(4>a — 4>b) 
2 cos 0a cos 9b)- 



To obtain the partition function, we further determine 
the integration measure. As two angles (9, 4>) specify a 
point onS2, integration measure can be taken as the sur- 
face area covering 52 x S2, such that 

Z = / cLOa / cLOb sin smOs 
Jo Jo 

I dcp A / d0 s exp(-/3V). (7) 



(i 



(i 



Next we will explain how we can obtain surface energy 
from this function. Our purpose is to obtain the inter- 
action energy between molecules that are placed across 
the virtual surface. Let us denote U and U' as molecules 
above that surface, and denote D and D' as molecules 
under the surface. This is shown in Fig. 2. Then we need 
to obtain the thermal expectation value of binding energy 
between two bulks divided by the virtual surface. 



<v> = n I f dMj2 y i 

U,D ^ ^ TT n 

1 n / d ^ u J d v 



U,D " " U,D 

0(£ v*W V uu> +E u,o Vun+E D^D' v dd> ) 



e ~0(Eu^u' v uu'+T.u,o V UD+En^O> V DD>) ] (g) 

where, dfxjj is the integration measure on the orien- 
tation of dipole moment of molecule U. Vxju 1 is the in- 
teraction energy between dipole moments of molecules U 
and U'. Other notations can be understood in the same 
way. We calculate the expectation value of the interac- 
tion energy between two dipoles, one is in U and another 
one is in D. Then we use the mean field approximation 
by rewriting unremarked other dipole moments to their 
mean values. Since there is no specific direction and we 
know that the water has no spontaneous electric polari- 
sation, we take these values to be zero. 



•V -. 



Virtual Surface 




t- 

Fig. 2 < Vud >'• dotted interactions are neglected. 
Molecules 



\ f B t 



(6) 



Interaction 

Fig. 3 nearest neighbour interactions. c?=n -1 / 3 : 
molecular distance. 



3 



Then we have 



with 



< V > ~ E / I dfi D VuD( 



-0Vui 



pair:U,D ' 

I J duu J dp D e~^ 

E < V UD > • 

pair:U,D 



(9) 



In a meanwhile we utilise Eq. only for surface 
molecules like in Fig. 3, and we calculate the surface en- 
ergy in this nearest neighbour approximation. Then the 
surface energy U sur is given by < Vab > as, 



_ <Vab> 

u sur — n 1 v i 



(10) 



where the negative sign is to make the surface energy 
positive. N is the surface density of molecules expressed 
by n. 

N=±= n 2 '\ 
d A 

Further < Vab > is expressed by the partition func- 
tion. 



Z = J dfi A J dp B exp[-f3V A B] , 



< Vab >-- 



dlogZ 
dp ■ 



(11) 



(12) 



Then we obtain 



TJ -l n 2/3 dl °9 Z 
v su.r — ~ 1 L 



(13) 



2" dp ' 
Our approximations are the followings. 

1. Rough mean field approximation. 

2. Only the nearest neighbour interaction. 

3. Neglecting entropy term for the surface tension 



III. CALCULATION OF PARTITION 
FUNCTION 

By calculating 

|>7T |>7T 

Z = / d$A \ d9s sin 9a sin 9 b 
Jo Jo 

x d(f> A d<t> B exp(-0V), (14) 
Jo Jo 

and by putting it into Eq. I|13|) . we obtain the surface 
energy in this rough approximation. First we carry out 
the integration. Let us define 



hip) 



2tt r 2ir 



., , d(j) A / d(j) B exp[-pcos(</>A - <Pb) 



(15) 



Pn 2 

p = z sin 9a sin#B. 



(16) 



We have 



h{p) 



i r 



where Jo is the 0-th Bessel function, and so Io{p) is 
the modified Bessel function. Io{p) can be expanded as, 



h{p) 



(P/2) 2 
^ (ml) 2 

m=0 y ' 



(18) 



From | p |< ^ ~ 3.4, this expansion converges in the 
normal temperature. And so the higher terms can be 
neglected approximately. By using this relation, we have 



Z = (2?r) 2 / d9 A / d9 B sin 9 A sin t 
Jo Jo 
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xJo( — — 3 sm0Asm.0B)e i7,e ° r .(19) 



Then we define 



a ee 2/3- 



(20) 



For r = d = r^ 1 / 3 , a(T = 273if) = 6.70, a(T = 
300K) = 6.09, a(T = 373K) = 4.90. By using the 
expansion of modified Bessel function, we obtain 



Z = (2n) 2 £ 



(q/4) 2 " 



d0 A / d6» B 



m =o (™ ! )- Jo Jo 
x sin 2m+1 6»Asm 2m+1 6l s e QCOS ^ coseB . (21) 

By setting x = con 9 a, y = cos$b, we obtain 



ro=0 v ' 

x(l-x 2 ) m (l-y 2 ) m e axy . 



We start from y-integration. 
r+i 



(22) 



dy(l - y 2 ) m e axy 



m / \ r+ l 

E (',') / x V) fc e- ' 



fc=0 



cr fc sinh(aa;) 
dx 2 * ax 



(23) 
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Then the partition function becomes 

oo 

Z = 2(2tt) 2 Z m (a), (24) 

m— 

And we look for the perturbative solution. 

2b (a) - P , (26) 



2 1 (a) = (^) 2 [P 1 --^P 2 ], 



(27) 



Z 2 (a) = ^) 4 [P 3 - 4,P 4 + -M (28) 



where, 



r+1 

P l = J dx { 1 ~ X 



ax a 
sinh(a.x) 



(29) 



ax 



2 2 2 
-Shi(a) H -sinh(a) ^cosh(a) (30) 



a J a 



dx{l - x 2 ) 



2 d sinh(aa;) 



eh; 2 



P* 



4 4 

--Shi(a) H sinh(a) 

a a 



2 2 2 

—Shi(a) ^(1 H ^)sinh(a) 

a a A a 

2 2 
r(l y)cosh(a), 



(31) 



(32) 



ci 2 sinh(ax) 
dx 2 ax 



P 4 = y da^l-s^)' 

8 24 24 
= Shi (a) ^sinh(a) H ~cosh(pt), (33) 

p 5 = f + \ x{ i- x ^^- s j^M 

J_i ax 4 ax 

48 64 

= — Shi(a) sinh(a) + 16cosh(a), (34) 

a a 



where Shi{a) is defined by 



. /" Q sinhz a 3 a 5 

Sfo(a) = y -— dz = a + — + — 



(35) 



Therefore, up to second order we have 

v 32a ' 32 ' 256 



Z (2) /(2tt) 2 = 2^Z„ = (^ + ^ + — )Shi(a) 



n=0 

+ ik )smh{a) - + h )cosh{a) - (36) 

Note that if we expand V(r) = —dlnZ/dp by a, we 
have 



< V(r) >iowest= ~ kT j28 C 



21 o 2 ^ 4 



3 fcT(47re) 2 r 6 ' 



(37) 



This expression of the lowest order in a is well known 
in many text books as weak field approximation. 
But in short range interaction (r ~ d), a expansion does 
not converge and such expression fails. In this sense our 
approximation(cxpansion of modified Bessel function) is 
essentially different from such a weak field approxima- 
tion. 

Now we obtain surface energy U sur and its temperature 
gvaaievAdU sur / dT up to second order in the following. 



where 



Z( 2 ) da 



dU, 



(2) 



dT 



-2k( 



47T60 71 



da 



(39) 



(40) 



(2) 

The surface energy Us Ur an d its temperature gradient 

( 21 

dUsur/dT at normal temperature 300A(a = 6.14) is 

U£l = 80.1[mJ/m 2 ], dU^ r /dT = -0. 108 [mJ/m 2 K]. 

(41) 

These values are comparable to the static values of the 
surface tension and its gradient, 7 exp = 73[mJ/m 2 ], 
{d 1 /dT) exp = -0.15[mJ/m 2 A]. 

Though this result looks good agreement with static 
value of surface tension, it is reported that, in experi- 
ments using oscillating jet method Q , the value of sur- 
face tension of water is about 180 [mJ/m 2 ] just after the 
creation of new surface @. And its value relaxes to the 
static value 73[mJ/m 2 ] within less than 10[mS]. This 
phenomenon is usually explained by the surface orienta- 
tion. The water molecules on the surface rotate and to 
make its surface free energy (surface tension) minimum 
during this period. Since we can not take into account the 
surface orientation from technical reason, we need to im- 
prove the approximation beyond the nearest-neighbour's 
one and to obtain the value 180[toJ/to 2 ] as the surface 
tension. This is discussed in the next section. 
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IV. BEYOND THE NEAREST-NEIGHBOUR 
APPROXIMATION 

As we have seen in the previous section, the form of 
binding energy is complicated even for one molecule-pair 
interaction, and it seems quite difficult to consider the 
theory beyond the nearest neighbour interactions. Fur- 
thermore, since the potential energy between two dipoles 
has the form of 1/r 3 , the sum of the potential energies 
diverges logarithmically. Therefore we give an important 
hypothesis, that is, we utilise thermal expectation value 
of dipole-dipole interaction energy as an elementary po- 
tential energy, and to sum up these energies for all the 
molecule pairs across the surface (See Fig. 5). 
From this hypothesis, the potential depends only on the 
distance of dipole pair. The function of potential starts 
from 1/r 6 , and so we have no divergence in summing up. 

From the previous calculation, we see that 
the partition function of dipole-dipole interaction 
can be expressed by even function of a. This 
means the expectation value of interaction energy: 
< V(r) >= —dlnZ/d/3 has the distance depen- 
dence in the series of l/r 6m (m = 1,2,3,...), since 
a = const, x f3/r 3 . 

Therefore we can write 

oc oo 

< V(r) >= kT B ma 2m = £ C m r~^ . (42) 

m—l m— 1 

For the explicit form of B m and C m , see appendix 1. By 
using above relation surface energy in nearest-neighbour 
approximation can be expressed as, 



U, 



N.N 



1 1 °° 



(43) 

Next we consider the case that all the molecules in a 
bulk interacts with all the ones in another bulk. 



Virtual Surface 



Bulk 



Interaction-. T distance L 
Molecule A 



Fig.4 
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Fig.5 



First as in Fig.4, the total energy that the molecule A 
apart L from the surface interacts with all the molecules 
in bulk is, 

/•oo />oo 

U(L) = n / dr2Trr / dz < V(y/r 2 + (z + L) 2 ) > 



E 



c„ 



- — 6m+3 



(3m — l)(6m — 3) 



(44) 



Second as in Fig.5, we sum up the energy U (L) by chang- 
ing the position of molecule A in opposite bulk. This 
summation is also calculated by using continuous approx- 
imation, i.e. by integration. The distance of two bulks 
e is an order of d, and its value is quite critical to the 
result. So we discuss its value later. Then we have the 
"bulk to bulk" surface energy. 



U 



B.B 



--n I dzU(e + z) 



9 OO 



C ... 



-6m+4 



^ (3m - l)(6m - 3) (6m - 4) 



(45) 



For e = d, this surface energy is smaller than the one 
in the case of nearest neighbour approximation, though 
it should be larger, (compare Eq. 143(1 with Eq. I|45() 
by noting n = d~ 3 .) The reason is the following. Since 
we have approximated all the molecules are continuously 
spread out, the bulk to bulk distance has increased effec- 
tively. This is shown in Fig. 6. Each molecule is spread 
out continuously in V± and Vi as shown in (b). For e = d, 
the mean distance between two molecules is larger than 
d. To determine the virtual bulk-bulk distance e, our 
consideration follows. 



o 



o 

(a) 



Vl l~ 

•tj 

v 2 [ --—- 



(b) 



Fig.6 



We require that the potential l/d 6m as in Fig.6 (a) co- 
incides with the mean potential < l/r 6m > as in Fig.6 
(b), where each molecule is spread out uniformly in col- 
umn V\ and V<z with hight d and basal area d 2 . Then we 
can determine e for each number of m. 



V. 



(m) 



I 

d^ 



VxV 2 



rf 3 



ri 



rf 3 



r-2 1 



V-2 



ri - r 2 



1 6m ' 



(46) 



Though it might seem strange that e depends on m, the 
bulk distance e is virtual mathematical object, and it 
improves the errors coming from continuous approxima- 
tion at short distance. Therefore it can be varied by 



6 



the parameter of potential form to. Hereafter we rewrite 
e — * e m- As we have shown in appendix 2, in large m 
limit this condition reduces to 



1 



7T 



1 



d 6m ~ 4 (3m- l)(6m-3)(6m-4) 



6m — 4 " 



(47) 



If we compare Eq. I(43|) with Eq. 1)45(1 by using above 
relation, we see the coincidence of the series expansion 
of U BB and U NN at large to. This is not surprising. 
For large to, the short range interaction enhanced, and 
the nearest neighbour interaction plays a central role in 
surface energy. For to = 1, 2 we obtain in appendix 2, 

ex = 0.373d, e 2 = 0.54639(1 

And for to > 3, Eq. I|47|l approximately holds and we 
just rewrite the series expansion of U BB by using U . 



U, 



N.N 



1 oo 



111 — 1 



TJ B - B 1 ^2^r < ^l ^-2 , -81 

^ - -2 Tl7r[ T2 £l + 360 62 ] 



1 



,2/3 



(48) 



m— 3 



We should remark here that these expansions are not 
convergent for real value of d. But we had to expand the 
potential energy in powers of r to sum up the interaction 
energies for all the molecule-pairs across the virtual sur- 
face. Though our two expansions are not convergent, the 
difference of these two surface energies is convergent. By 
using the value of e which are calculated in appendix 2, 
we have 



U. 



B.B 



N.N 

sur 



jjN. 



l n2/ 3 [ Q^ + 0.09856a 2] (49) 



The r.h.s. is the correction to nearest-neighbour approx- 
imation, which comes from distant interactions. From 
C m /d 6m = kTB m a 2m , we have 



U, 



B.B 



u. 



N.N 



2(3 



[0.8817Sia 2 + 0.09856B 2 a 4 



4ire 



[0.8817Sia + 0.0985652a 3 
21 



= - 131 [0.8817 x(-— )a 
9779 

+0 - 09856 x (2457600^ 3 ]^- 



(50) 



Since we know the analytical value of U^J r , all the 
surface energy due to dipole-dipole interaction becomes 

V™ = - (0-8817 x (-|L )a 

9779 

+0 ' 09856 X (2457600^ 

= 131 x [/(a) + 0.14465a - 0.0003922a 3 ], (51) 



where the unit is given by [mJ]. The Fig. 7 shows the a 
dependence of the surface energy. 



Usur [mj/m' 




a 



Fig. 7 U s . 



a 



Our theory gives the surface energy 185[toJ/to 2 ] at 
normal temperature (300i"Q. This result well coincides 
with the experimental value 180[toJ/to 2 ] by using the 
oscillating jet method @, 0- However, we need to take 
into account other interactions. This is done in next sec- 
tion. 



V. ADDITIONAL OTHER FORCES 

We still have 3 kind of forces that work between 
molecules such as, dispersion force, induced force, and 
repulsive force (part of Lennard- Jones potential). So we 
must include them into our consideration of surface en- 
ergy. 

The dispersion force is essentially quantum mechanical 
effect. 



TT - 3 ^ ft ° ^ J 
dlsp ~ 4W j r 6 ' 



(52) 



where ao, the electronic polarisation rate and /, the 
first ionization energy. These values are given as 
a /(47re ) = 1.48 x 10 - 30 [to 3 ], / = 12.6[eF]. 

The induced force is the interaction between static 
molecular dipole and the instantaneous dipole moment 
due to the fluctuation of electron cloud. [4| 



Uind — 



2M 2 qp 

"(47re ) 2 r 6 ' 



(53) 



The experimental value of the repulsive energy, a part 



of Lennard-Jones potential is, 



le repuJ 

MM 



U re p 4s ( ) , 
r 

where e/k = 809.1[if], a = 2.641 x 10~ 10 [to]- 



(54) 



We sum up these 3 energies by using d = n ^ 
3.224 x 10- 10 [to], 



AC/ = U dl 



sp 



U, 



ijid 



u, 



■rep 



-4.97(-) 6 + 4.08(-) 12 ] x 10~ 18 [toJ]. (55) 
r r 
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Now our consideration goes to the surface energy as be- 
fore. We sum up all the interaction energies across the 
virtual interface. For the molecule-molecule potential 



Bulk-bulk surface energy becomes 



U sur — ~ n> ^ ^ 



-6m+4 



^ (3m - l)(6m - 3)(6ra - 4) 



where t\ = 0.373d, e 2 = 0.5463M. In [mJ] unit we 
obtain for this additional energy 



Ci = -4.97 x 10 



-18 j6 



C-2 = 4.08 x lO -18 ^ 



Then we have additional surface energy due to these 3 
forces. 



4.97, 1 



4 .08 . 



12 v 0.373 



1 



360 '0.54639 



) 8 ] x 10~ 18 = 25[toJ/to 2 ]. (56) 



From above calculations, the total surface energy at 
normal temperature is 

U sur = U dipole + U dtsp + U md + U rep = 210. (57) 

Comparing to the experimental value of surface tension: 
180[mJ/m 2 ], the obtained value is only 17 % larger. [|| 
So we may conclude that our method to obtain the sur- 
face tension expresses the nature of water. 



In the same way, 



< V(r) > _^ ^<91ogZ (2) 



kT 



00 



21 2 9779 4 
= -T^a + — or 



128 2457600 
650953 



92250213859 



4404019200 15981304872960000 
7557351481321 



33002459982987264000 



(59) 



Each coefficient of a 2m is B m . Then the form of C m is 
given by 



Cm — P B m {— ) 

47re 



VIII. APPENDIX 2 



Vl 



Interaction 



V2 




(60) 



VI. SUMMARY 



Fig- 



We have calculated the dipole-dipole interaction en- 
ergy due to new approximation by using statistical me- 
chanics. Then we sum up the energies for all the molecule 
pairs across the surface in bulk. In this way we have cal- 
culated the surface energy of water as dipoler liquid. Our 
treatment is the first step for obtaining the surface en- 
ergy, though some improvements may be necessary. The 
first improvement should be done for the rough mean field 
approximation. The thermal expectation value < Vjjd > 
may be affected by other environmental molecules. This 
might be done by introducing the effect of dielectric con- 
stant like the Lifschitz theory. 



VII. APPENDIX 1 



We expand the partition function Eq. (|36|l by a series. 



21 



91 



%)/(^ 2 =4 + ^ 2 + 960(r 



113 

"806400* 



407 



304819200 



a° + .. 



(58) 



y(r. 



1 



dL 



1 

7p 



dz 



2irrdr 



y/(e + L + z) 2 + r 



-6m " 



(61) 



Since we have a translational invariance to the inter 
facial direction, when we integrate by r\ in column V\ 1 
we fix ?i in the centre axis in column Vi- After that, the 
integration of ?i is given only in the centre axis direction 
such as, integrating from to d and divided by d. 

First r-integration can be done and we obtain 



y( m ) = 



(3m 



1 



[d 2 /TT+{e + L + z) 2 ] 3m 



(e + L + zf m - 2 
IT)- (62) 
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Next by using the formula 



dx I dyg(x + y) 
'a Jo 

— 2 / du g(u) + / du u g(u) — I du u g(u), (63) 
Ji Jo Ji 

and by setting x = z/d, y = L/d, u = x + y, e' = e/d, 
y( m ) = 



(3m- l)tf 



(ill! 



,2+e' 1 1 

Jl+e' ^tfl^ ~ (l/TT + U 2 ) 3 ™- 1 ^ 

2+£ ' , 1 1 

1+6 , ^ _ £ ')(^^2 - ( 1/7r + M 2)3rn-1 ) 

1+e' 



+ / du(u — e')( 



1 



u 6m-2 (I/tt + m 2 ) 3 " 1 - 



-)]-(64) 



In the first and second integration term, the integrand 
becomes smaller since u > 1 for large m, Therefore we 
have 



y(m) I 



(3m - l)d 6m 
1 1 



1+e' 



du(u — e ) 



E [ I 

(3m- l)d 6mL (6m-4)(6m-3)e' 6m - 4 

'' 1+£, u - e' 

C?M- ; ^T"5 r l. (65) 

(l/7T + U 2 ) 3m - lJ V ^ 

From the numerical calculation, the rate of the last inte- 
gration term to the first term is about 8% for m = 2 and 
2.5% for m = 3. So we can neglect the last integration 
term for m > 3. Therefore the condition l/d 6m = Ve 
is simplified in large m limit as 



(3m-l)(6m-3)(6m-4) 1 



(66) 



If we compare Eq. (|43() with Eq. I|45(l using above re- 
lation, we see the coincidence of the series expansion of 
U BB and JJ NN at large m. This is quite natural result. 
Because at large m, only the short range interaction en- 
hanced, and so the nearest neighbour interaction gives 
the main effect. 



The condition l/d 6m = Ve" 1 '* is now rewritten 



(3m-l)(6m-3)(6m-4) 1 



as 



7r e 
1 2 



/6m— 4 



(2 + e ')6m-4 ^ + £ /)6m-4 

-(6m — 3)[ 



((2 + e') 2 + l/7r) 3 '"- 2 
2 1 



((l + e') 2 + iA) 3m ~ 2 (e' 2 + l/7r) 2 



- ( 6m-3) ( 6m-4)[ ( 2 + e0y i+e/ ^ + 



2+e' 



1 



du 



1 



(w 2 + l/7r) 3m - ir 



(67) 



The remained integrations in right hand side can be cal- 
culated in series expansion for general m. For m = 1, 
this condition becomes 



12 1 1 2 

~ ~ ~7> ~ (2 + e') 2 ~ (1 + e') 2 

+3V^ 3 [-(2 + e') arctan(V^(2 + e')) 

+2(1 + e') arctan(V^(l + e')) - e' arctan(V7i : e / )] (68) 



By solving this, we have = 0.373. In the same way we 
obtain e' 2 = 0.54639, and e' 3 = 0.63718926. If we use Eq. 
(|66|l for m = 3 instead, we have £3 = 0.638. From above 
calculation we use above numerical value for m = 1,2, 
and we use U NN instead of U BB for m > 3. 
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